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Approximation for non-smooth functionals of stochastic 
differential equations with irregular drift 


Hoang-Long Ngo* and Dai Taguchr 


Abstract 

This paper aims at developing a systematic study for the weak rate of convergence of 
the Euler-Maruyama scheme for stochastic differential equations with very irregular drift and 
constant diffusion coefficients. We apply our method to obtain the rates of approximation 
for the expectation of various non-smooth functionals of both stochastic differential equations 
and killed diffusion. We also apply our method to the study of the weak approximation of 
reflected stochastic differential equations whose drift is Holder continuous. 
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1 Introduction 

Let T > 0 be fixed and (X t )o <t<T be the solution to 

dX t = b(X t )dt + <r(X t )dW t , X 0 = x 0 £R d , 0 <t<T, 

where W is a d-dimensional Brownian motion. The diffusion ( X t )o<t<T is used to model many 
random dynamical phenomena in many fields of applications. In practice, one often encounters the 
problem of evaluating functionals of the type E[/(V)] for some given function / : C[0, T} —> M. For 
example, in mathematical finance the function / is commonly referred as a payoff function. Since 
they are rarely analytically tractable, these expectations are usually approximated using numerical 
schemes. One of the most popular approximation methods is the Monte Carlo Euler-Maruyama 
method which consists of two steps: 

1. The diffusion process (X t )o<t<T is approximated using the Euler-Maruyama scheme (X d ) 0 < t <T 
with a small time step h > 0: 

dXt = b{X* h{t) )dt + a(X* h(t) )dW u X% = x 0 , rj h (t) = kh, (1) 

for t £ [kh, ( k + 1 )h), k £ N. 
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2. The expectation E[/(X)] is approximated using f(X h ’ x ) where (X h ’ x )i<i <jv are N 

independent copies of X h . 

This approximation procedure is influenced by two sources of errors: a discretization error and a 
statistical error 

1 N 

Err{f , h) := £rr(fc) := E[/pQ] - E [f(X h )}, and E[/(X h )] - - £ 

V 2=1 

We say that the Euler-Maruyama approximation ( X h ) is of weak order re > 0 for a class 7~L of 
functions / if there exists a constant K(T ) such that for any f G 71, 

\Err(f,h)\<K(T)h K . 

The effect of the statistical error can be handled by the classical central limit theorem or large 
deviation theory. Roughly speaking, if f(Xlf) has a bounded variance, the L 2 -norm of the statistical 
error is bounded by N~ 1 ^ 2 Var(X^) 1 ^ 2 . Hence, if the Euler-Maruyama approximation is of weak 
order re, the optimal choice of the number of Monte Carlo iterations should be N = 0(h ~ 2K ) 
in order to minimize the computational cost. Therefore, it is of both theoretical and practical 
importance to understand the weak order of the Euler-Maruyama approximation. 

It has been shown that under sufficient regularity on the coefficients b and a as well as /, 
the weak order of the Euler-Maruyama approximation is 1. This fact is proven by writing the 
discretization error Err(f , h) as a sum of terms involving the solution of a parabolic partial dif¬ 
ferential equation (see GO ES m H3 Q3])- It should be noted here that besides the Monte Carlo 
Euler-Maruyama method, there are many other related approximation schemes for E[/(Xy)] which 
have either higher weak order or lower computational cost. For example, one can use Romberg 
extrapolation technique to obtain very high weak order as long as Err(h) can be expanded in terms 
of powers of h (see (S3). When / is a Lipschitz function and the strong rate of approximation 
is known, one can implement a Multi-level Monte Carlo simulation which can significantly reduce 
the computation cost of approximating E[/(X)] in many cases (see [9]). It is also worth looking at 
some algebraic schemes introduced in [25j. However, all the accelerated schemes mentioned above 
require sufficient regularity conditions on the coefficients b,a and the test function /. 

The stochastic differential equations with non-smooth drift appear in many applications, es¬ 
pecially when one wants to model sudden changes in the trend of a certain random dynamical 
phenomenon (see e.g., [20]). There are many papers studying the Euler-Maruyama approxima¬ 
tions in this context. In [JO] (see also [3]), it is shown that when the drift is only measurable, 
the diffusion coefficient is non-degenerate and Lipschitz continuous then the Euler-Maruyama ap¬ 
proximations converges to the solution of stochastic differential equation. The weak order of the 
Euler-Maruyama scheme when both coefficients b and a as well as payoff functions / are Holder 
continuous has been studied in [201 EL]. In the papers |2T[ and j33j [34] , the authors studied the 
weak and strong convergent rates of the Euler-Maruyama scheme for specific classes of stochastic 
differential equations with discontinuous drift. We would like to note here some very recent papers 
on the unbiased simulation [T6] and exact simulation in 0 0IM] methods for stochastic differential 
equations with irregular coefficients. 

The aim of the present paper is to investigate the weak order of the Euler-Maruyama approx¬ 
imation for stochastic differential equations whose diffusion coefficient <r is constant, whereas the 
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drift coefficient b may have a very low regularity, or could even be discontinuous. More precisely, 
we consider a class A of functions which contains not only smooth functions but also some dis¬ 
continuous one such as indicator function. The drift b will then be assumed to be either in A or 
a-Holder continuous. It should be noted that no smoothness assumption on the payoff function 
/ is needed in our framework. As a by product of our method, we establish the weak order of 
the Euler-Maruyama approximation for some particular functionals / which include the path-wise 
maximum of the diffusion, integral of diffusion with respect to time as well as the approximation 
of a diffusion processes killed when it leaves an open set. We also apply our method to study the 
weak approximation of reflected stochastic differential equation whose drift is Holder continuous. 

The remainder of this paper is organized as follows. In the next section we introduce some 
notations and assumptions for our framework together with the main results. All proofs are 
deferred to Section [3] 


2 Main Results 


2.1 Notations 


For an invertible d x d-matrix A = (Aij)i<ij<d, we define for all x,y £ R d , 


9A(x,y) 


exp(-^(A 1 (y - x),y - x)) 
(27r) d / 2 v / det ~A 


( 2 ) 


In particular we denote g c (x, y) = g c i{x, y ) for c £ R where the matrix / is the identity matrix. 

A function £ : R d —> R is called exponentially bounded or polynomially bounded if there exist 
positive constants K,p such that |£(a;)| < Ke K \ x I or |£(x)| < K( 1 + |x| p ), respectively. 

Let A be a class of exponentially bounded functions £ : R d —> R such that there exists a 
sequence of functions ((n)ngn C C 1 { R d ) satisfying: 


MW: Gv -> £ in Lj oc (R d ), 

< A{ii) : supjv |£jv(a:)| + |£(x)| < Ke K ^, for all x G R d (3) 

[A{iii) : supjv u>0; oeR de- K M- Ku f Rd |V£jv (x + a)\dx < K, 

for some positive constant K. We call (£/v)jveN a,n approximation sequence of £ in A. 

The following propositions shows that this class is quite large. 

Proposition 2.1. i) If £, £ € A then ££ £ A and ai£ + a 2 £ € A for any ai,a 2 € M. 

ii) Suppose that A is a non-singular d x d-matrix, B £ R d . Then £ £ A iff f(x) := £(Ax + H) £ A, 

for all x £ R d . 


It is easy to verify that the class A contains all C' 1 (R d ) functions which has all first order 
derivatives polynomially bounded. Furthermore, the class A contains also some non-smooth func¬ 
tions of the type £(a;) = {x\ — a) + or £(x) = I a <x<b for some a,b £ R d . Moreover, we call 
a function £ : —> R monotone in each variable separately if for each i = 1 ,...,d, the map 

Xi i y £(xi,..., Xi,..., x n ) is monotone for all xi ,..., Xi-\,Xi+i ,..., Xd £ R. 


3 





Proposition 2.2. Class A contains all exponentially bounded functions which are monotone in 
each variable separately. 

The proofs of Propositions 12.11 and 12.21 and further properties of class A were presented in YU1\ 
and [531 . 

We recall that a function £ : R d —>- R d is called a-Holder continuous for some a £ (0,1] if 
there exists a positive constant C such that |C(a;) — C(d)l < C \x — y\ a for all x, y £ R d . We denote 
A4(R d ,R d ) the class of all measurable functions b : R d —> R d . Let a £ (0,1] and let B(a) be the 
class fo functions given by 

B(a) = {be 7W(R d ,R d ); b = b H + b A where b H is a-Holder and bf £ A, j = 1,..., d) . (4) 

A function b £ B(a) could be very irregular since its trajectory may both change continuously 
very fast because of the Holder part b H , and be discontinuous at some points because of the part 
b A £ A. Moreover, the class A does not contain any class of Holder continuous functions H 13 for 
any /3 £ (0,1) (see j5il). 


2.2 Weak approximation of stochastic differential equations 


Let (B, X, (-7 r t)t>0;I I> ) be a filtered probability space and (W t ) t >o be a d-dimensional standard 
Brownian motion. We consider a d-dimensional stochastic differential equation 

X t = x 0 + f b{X s )ds + aW t , Xq £ R d , t £ [0,T], (5) 

do 

where a is a d x d deterministic, uniformly elliptic matrix, that is for the matrix a := era*, there 
exist 0 < a < a < oo such that for any £ £ R d , 


and b : — > R d is a Borel measurable function. Let X h ,h > 0, denote the Euler-Maruyama 

approximation of X , 




xo + [ b(XK s) )ds + aW t , t £ [0,T], 


( 6 ) 


where r)h(s) = kh if kh < s < (k + 1 )h for some nonnegative integer k. In this paper, we study the 
convergent rates of the error 

Err(h) = E[f(X)] - E[f(X h )} 
as h —> 0 for some payoff function / : (7[0, T] —> R. 


Remark 2.3. Note that the uniformly elliptic condition plays an important role in establishing the 
convergence of the Euler-Maruyama approximation for SDEs with non-Lipschitz coefficients. In 
fact, Hairer, Hutzenthaler and Jentzen m Theorem 5.1] constructed a class of A-dimensional SDEs 
whose drift coefficient is a smooth, bounded and non-Lipschitz function, and diffusion coefficient 
is a deterministic non-uniformly elliptic matrix for which the Euler-Maruyama scheme does not 
converges with any polynomial rate, that is 


lim E[\X T -X*\] = lim \E[X T ]-E[X*]\ = 
h \o h a h\0 h a 
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0 

oo 


ifa = 0, 
if a > 0. 




A Borel measurable function £ : R d —> R d is called sub-linear growth if £ is bounded on 
compact sets and £(y) = o(|y|) as y —> oo. £ is called linear growth if |£(y)| < c\\y\ + C 2 for some 
positive constants cr,C 2 - 

Remark 2.4. (i) It has been shown recently in that when b is of super-linear growth, i.e., 

there exist constants C > 0 and 9 > 1 such that |£(y)| > \y\ e for all \y\ > C, then the 
Euler-Maruyama approximation © converges neither in the strong mean square sense nor 
in weak sense to the exact solution at a finite time point. It means that */E[|Xt’| p ] < oo for 
some p 6 [1, oo) then 

lim E \\X T - X£| p l = oo and lim |E[|X t | p - |X£| P 11 = oo. 
h-s-O L J 1 L J 1 

Thus, in this paper we will consider the case that b is of at most linear growth. 

(ii) In the one-dimensional case, d = 1, it is well-known that if a ^ 0 and b is of linear growth, 
then the strong existence and path-wise uniqueness hold for the equation w (see W- 

(in) In the multidimensional case, d > 1, it has been shown in m that if b is bounded then the 
equation 0) has a strong solution and the solution of 0) is strongly unique. Moreover, if 
a is the identity matrix, then the equation 0 has a unique strong solution in the class of 
continuous processes such that P (f 0 |6(X s )| 2 ds < oo) = 1 provided that f Rd \b(y)\ p dy < oo 
for some p > d V 2 (see { 2 ^). 

Throughout this paper, we suppose that equation (0 has a weak solution which is unique in 
the sense of probability law (see Chapter 5 jl9|). 

The following results requires no assumption on the smoothness of /. 

Theorem 2.5. Let X be a solution of the SDE 0. Suppose that b G £>(a) and b is of linear 
growth. Moreover, assume that f : C[0,T] —> R is bounded. Then 

hmE[f(X h ))=E[f(X)]. 

h— s-0 


If b is of sub-linear growth, we can obtain the rate of weak convergence as follows. 

Theorem 2.6. Let X be a solution of the SDE 0 and h > 0. Suppose that b € B(a) and b is 
of sub-linear growth. Moreover, assume that f : C[0,T] — > R satisfies E[|/(xo + crIT)| r ] < oo for 
some r > 2. Then there exists a constant C = C(xo,b,a,T), which depends neither on h nor on 
f, such that 

|E lf(X)} - E[f(X h )}\ < c(E[\f(x 0 + aW)\ r ]) Vr h^i. 

The following result concerns with the approximation of maximum of SDEs. 

Corollary 2.7. Assume the hypotheses of Theorem \ 2. 61 Moreover, suppose that g : R —> R is 
/ 3-Holder continuous with /3 € (0,1]. Then there exists a constant C which does not depend of h 
such that 


E 


9 


( max I X s I 

0 <s<T 


-E 


9 


max 

V°< S < T 



<C 


(hlog(l/h))^ 2 } . 
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For an integral type functional, we obtain the following corollary. 

Corollary 2.8. Let h = T/n for some n G N. If the drift coefficient b G 13(a) is bounded, then for 
any Lipschitz continuous function f and g G B(/3) with ft G (0,1], there exists a constant C which 
does not depend of h such that 


E 


g(X s )ds 



< C7j* A * A *. 


Remark 2.9. In the paper !2df . the author considered the weak rate of convergence of the Euler- 
Maruyama scheme for equation w in the case of a one-dimensional diffusion. It was claimed that 
if b was Lipschitz continuous, the weak rate of approximation is of order 1. However, we would 
like to point out that the given proof contains several gaps (see for instance Lemma 2 of J29/ and 
Remark \3.3\ below 1 which leave us unsure about the claim. 

Remark 2.10. It has been shown in [H]/ that for a stochastic differential equation with a- 
Holder continuous drift and diffusion coefficients with a G (0,1), one has 

\E[f(X T )]-E[f(X!f)]\<Ch a / 2 , 

where f G C 2 and the second derivative of f is a-Holder continuous. On the other hand, in m 
Gyongy and Rasonyi have obtained the strong rate of convergence for a one-dimensional stochastic 
differential equation whose drift is the sum of a Lipschitz continuous and a monotone decreasing 
Holder continuous function, and its diffusion coefficient is Holder continuous. In 1 33/ . the authors 
improve the results in m More precisely, we assume that the drift coefficient b is a bounded and 
one-sided Lipschitz function, i.e., there exists a positive constant L such that for any x,y G R d , 
(x — y , b(x) — b(y)) K d < L\x — y\ 2 , bj G A for any j = 1,..., d and the diffusion coefficient a is 
bounded, uniformly elliptic and (1/2 + a)-H6lder continuous with a G [0,1/2]. Then for h = T/n, 
it holds that 


( C(logl//i) 1 if a = 0 and d = 1, 

E\\X T - X!f\] < I Ch a if a € (0,1/2] and d= l, 

\ Ch 1/2 if oi= 1/2 andd> 2. 

Therefore, if the payoff function f is Lipschitz continuous, it is straightforward to verify that 

( C(log l/hff 1 ifa = 0andd = l, 

|E [f(X T ) - f(X$)\ | <1 Ch a ifae (0,1/2] and d= 1, 

[ Ch 1 / 2 if a = 1/2 and d> 2. 

In the following we consider a special case of the functional /. More precisely, we are interested 
in the law at time T of the diffusion X killed when it leaves an open set. Let D be an open subset 
of R d and denote td = inf(f > 0 : X t £ D}. Quantities of the type E[g(XT)l( TD >T)\ appear in 
many domains, e.g. in financial mathematics when one computes the price of a barrier option on 
a d-dimensional asset price random variable X t with characteristics /, T and D (see mm and 
the references therein for more detail). We approximate td by = inf {kh > 0 : X( l h £ D,k = 
0 , 1 ,...}. 
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Corollary 2.11. Assume the hypotheses of Theorem \2.6\ Furthermore, we assume 
(i) D is of class C°° and dD is a compact set (see fffj and IlOf ): 

(ii) g : R d —> R is a measurable function, satisfying d(Supp(y), dD) > 2e for some e > 0 and 
llfflloo = SUp^gRd |ff(x)| < OO. 

Then for any p > 1, there exist constants C and C p independent of h such that 


E[ff(Xr) 1 (t d >T)\ E[a(^T')l(r^ > T)] 


< Ch 2A I 


Cpllg p lloo 

1 A e 4 /p 


h 2 ^. 


(7) 


Finally, we consider the approximation for the density of SDE ([5]). Let pt(x o, •) and p/(x o, •) 
be the density functions of X t and Xjf respectively. Then we have the following rate of convergence. 

Theorem 2.12. Suppose that b £ 13(a) and bounded. Then for any p > d and r > 1, there exists 
constants C PtT and c p such that for any y £ and h £ (0,T/2), 

\pr(x 0 ,y) - Pr(xo,y)\ < C p , r gc p T(x 0 ,y) + h a/2 + h 1/{2pr) j , 

where g Cp T(xo,y) is the density function of normal distribution defined in ©. 

Remark 2.13. (i) Note that if d = 1, we can choose p £ (1,2) and r = 2/p, and then 

\pt(x 0 , y) - pt(xo , y) I < C Ptr g Cp T(x 0 , y) | + /jl/4 } • 

(ii) Konakov and Menozzi 1 2Sj obtain a better rate of convergence under a further assumption that 
the drift coefficient is piecewise smooth (see Theorem 2 in 1 23/ ). However, in our setting, the drift 
coefficient may have infinite number of discontinuous points. 


2.3 Weak approximation of reflected stochastic differential equations 

We first recall the Skorohod problem. 

Lemma 2.14 ([T9], Lemma III.6.14). Let z > 0 be a given number and y : [0, oo) — > R be a 
continuous function with yo = 0. Then there exists a unique continuous function t = (lt)t>o 
satisfying the following conditions: 

(i) x t •■= z + y t + £ t > 0, 0 < t < oo; 

(ii) i is a non-decreasing function with = 0 and it = /q 1(x s = 0 )di s . 

Moreover, i = (it)t >o is given by 

it = max{0, max (—z — y s )\ = max max(0, i s — x s ). 

0<s<t o <s<t 


7 







Let us consider the following one-dimensional reflected stochastic differential equation valued 
in [0, oo) such that 


(x t — xo + f Q b(X s )ds + aWt + L^{X), Xq £ [0, oo), t £ [0, T], 

|L°(X) =f*l (Xs=0) dL°(X), 

where (Lf(X))o<t<r is a non-decreasing continuous process stating at the origin and L°(X) is 
called the symmetric local time of X up to time t at the level 0. Equation © will be studied 
under the assumptions that b is of sub-linear growth and Holder continuous and a is a non-zero 
constant. Moreover, we assume that the SDE © has a weak solution and the uniqueness in the 
sense of probability law holds (see j3H 00]) . Using Lemma 12.141 we have 

L°(X) = sup max (0, L° S (X) — Xf) . 

o <s<i 

Now we define the Euler-Maruyama scheme X h = (Xjf)o<t<T for the reflected stochastic 
differential equation ©. Let X[} := x 0 and define 

X t h =*o + f b(X% h(s) )ds + aW t + L° t (X h ). (9) 

J o 

The existence of the pair (Xf,Lf(X h ))o<t<T is deduced from Lemma 12.141 Moreover 

L° t (X h )= /‘l w=0 )dL° s (X h ). 

Jo 

By the definition of the Euler-Maruyama scheme, we have the following representation. For each 

k = 0,l ,..., 


X^k+i)h ^ Xk h + b(X% h )h + cr(W( k+ i) h — Wkh) + max(0, A k — X kh ), 


where 


A k := sup (~b(X£ h )(s - kh) - a(W s - W kh j) . 

kh<.s<(k-\-l)h 

Though A k is defined by the supremum of a stochastic process, it can be simulated by using the 
following lemma. 

Lemma 2.15 ([27], Theorem 1). Let t £ [0,T] and a, c £ R. Define S t := sup 0<s<t (aW' s + cs). 
Let Ut be a centered Gaussian random variable with variance t and let Vt be an exponential random 
variable with parameter 1/(2 1) independent from U t - Define 

Yt := —(aUt + ct + (a 2 V) + (aUt + cf) 2 ) 1 / 2 ). 

Then the vectors ( Wt,St ) and ( Ut,Y t ) have the same law. 






Under the Lipschitz condition for the coefficients of the reflected SDE (JSJ), Lepingle j2S] shows 

that 


E[ sup \X t - X t h I 2 ] 1 / 2 < Ch 1/2 , 

0 <t<T 


for some constant C. 

We obtain the following result on the weak convergence for the Euler-Maruyama scheme for 
a reflected SDE with non-Lipschitz coefficient. 

Theorem 2.16. Let X be a solution of the SDE 0 and h > 0. Suppose that the drift coefficient 
b is of sub-linear growth and a-Holder continuous with a G (0,1]. Moreover, assume that f : 
C[0,T] — > R is bounded. Then there exists a constant C, which does not depend on h, such that 

\E[f(X)\-E[f(X h )]\<Ch a / 2 . 


3 Proofs 

From now on, we will repeatedly use without mentioning the following elementary estimate 

sup |afe fc|x| - x2 < oo, for any p > 0, k£ R. (10) 

Throughout this section, a symbol C stands for a positive generic constant independent of the 
discretization parameter h, which nonetheless may depend on time T, coefficients b , a and payoff 
function /. 

3.1 Change of Measures 

Let X be a solution of the SDE (0 defined in a certain probability space (12, J-, P) with a Brownian 
motion (Wt,Tt)t> o- The notations xq, b, er, rjh are defined as at the beginning of Section [2721 From 
now on, we will use the following notations 

Z t =e Yt , Y t = f (a~ 1 b) j (x 0 + crW s )dWf - ^ [ |cr _1 6(xo + crW s )| 2 ds, 

Jo * Jo 

Z h t = e Y ", Y t h = j\a~ 1 b) j (xo + aW^dWf - i j* \a~ l b{x 0 + <jW Vh(s) )\ 2 ds , 

where we use Einstein’s summation convention on repeated indices. We also use the following 
auxiliary stopping times 

rjf = inf{t > 0 : xo + crWt £ D}, and r^’ h = inf {kh > 0 : xo + crWkh & D, k = 0,1,...}. 

Lemma 3.1. Let X be a solution of the SDE 0 and h > 0. Suppose that b is a function with at 
most linear growth, then we have the following representations 

E[f(X)} - E [f{X h )\ = E[f(x 0 + aW)(Z T - Ztf% (11) 
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and 


= ®%(£ 0 + aW T ){Z T l( T w >T) - Z^l (r w,h >T) )], (12) 

for all measurable functions f : CfO, T] —> R and g : R d —> R provided that all the above expectations 
are integrable. 

Proof. Let <j~ 1 be the inverse matrix of er. Since b is of linear growth, so is a~ 1 b. Thus, there exist 
constants ci,C 2 > 0 such that |6(x)| < Ci|x| + C 2 for any x £ R d . For any 0 <t < to <T, 

\Xt\ < |a; 0 | + |crfFt| + f |6(X s )|ds 

Jo 

<|x 0 | + |cr| sup \W S \ + c 2 t 0 +Ci / |A s |ds. 

0<s<to JO 

Applying Gronwall’s inequality for t, £ [0, io]* one obtains 

|A to | < (|x 0 | + \cr\ sup \W S \ + c 2 t 0 )e Clt ° 

0<s<to 

< {\x 0 \ + c 2 T)e ClT + \<j\e ClT sup \W S \. (13) 

0<s<to 


On the other hand, for each integer k > 1, one has 


\Xkh\ < j + mxtk-DH )\+ m 


sup \W t \ 

0 <t<kh 


< (1 + hci)\Xk_ 1)h \ + hc 2 + 2|cr| sup \W t \. 

0 <t<kh 


Hence, a simple iteration yields that 


\^kh\ - (1 + hci) k \xo\ + (hc 2 + 2|cr| 


sup 

0<t<kh 


\w t \) 


(l + feci)^ 1 
hc\ 



Thus, for any t £ (0, T], 


!*£„(*) I < (l + /i Cl ) T/,l |a;o| 


c 2 (l + hci) T/h (1 + hci) T / h 


Cl 


+ 2H 


hc\ 


sup \W„ 

0<s<rj h (t) 


Moreover, 

\X^-X^ h{t) \< Cl h\Xl {t) \+c 2 h + 2\a\ sup \W t \. 

0<s<t 

Therefore, for any t £ (0, T], we have 


\X?\ < (l + Ci/i) 1+ ^ Cl| "° l + C2 +c 2 h + 

Cl 


2|cr| (1 + hci) 1+T / h + 2hc\ 
hci 


sup \W S \. 

0 <s<t 


(14) 
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We define new measures Q and Q h as 


dQ 

d¥ 


dQ h 


dP 


= 6XP ( ~~ Jo ^ lb ') j ( X ^ dW * ~ \ 

■ = exp ( - j\<r- 1 b) j (X* hlt) )dWi 


| a 1 b(X s )\ 2 ds'j, 

~\j Q I a ~ lb ( X v h (s))\ 2 ds)- 


It follows from Corollary 3.5.16 in [19] together with estimates m and m that Q and Q h 
are probability measures. Furthermore, it follows from Girsanov theorem that processes B = 


m 


B “), 0 < t < T} and B h = {(B 


h, 1 


Bt’ d ), 0 < t < T} defined by 


B[ = Wj + f*(a~ 1 b)j(X s )ds, B™ = Wj + 

Jo Jo 


1 < j < d, 0 < t < T, 


are d-dimensional Brownian motions with respect to Q and Q h , respectively. Note that X s = 
xq + crB s and X£ = xq + aB £. Therefore, 


E[f(X)] = Eq f(X) 


dP 


= Err 


/{xq + aB) exp ^ J (a 1 b) j (x 0 + crB s )dB J s - i J \a 1 b(x 0 + crB s )\ 2 dsj 


= E [f(x 0 + aW)Z T ]. 


Repeating the previous argument leads to E[/(X ?l )] = E[/(xo + aW)Zj], which implies (fTTTl . The 
proof of (fT^ll is similar and will be omitted. □ 


From now on, we will use the representation formulas in Lemma 13. II to analyze the weak rate 
of convergence. We need the following estimates on the moments of Z and Z h . 

Lemma 3.2. Suppose that b is of sub-linear growth. Then for any p > 0, 

E[\Z t \ p + \Ztf\ p ] < C < oo, 

for some constant C which is not depend on h. 


Proof. It suffices to proof the statement for p > 1. Using Holder’s inequality, we have 

r T r T 

E[e plT ] = E exp ^ p J (a~ 1 b)j(xo + aW s )dW° s — ^ J \a~ 1 b(x + crW,j)| 2 ds^ 

exp (p f (a^ 1 b)j(x 0 + <rW s )dWf - p 2 f \a~ 1 b(x 0 + crH4)| 2 ds+ 

v Jo Jt„-i 


= E 


+ ( V 2 - |) / W 1 b(x o + o-W s )| 2 ds) 


r T r T 

< |e exp ^2 p J (a~ 1 b)j(xo + aW s )dWj ^ 2p 2 J |cr _1 6(xo + trW s )| 2 ds^ J | 


1/2 


11 












pT 

|e exp((2p 2 -p) J |cr _1 6(x 0 + aW s )\ 2 ds^j j 


1/2 


Since b is of linear growth, so is er 1 b and it follows from Corollary 3.5.16 in [193 that 


exp ^2 p J ( a 1 b) J (x 0 + aW s )dW 3 s — 2p 2 J \a 1 b(xo + oW s )\ 2 dsJ =1. (15) 

On the other hand, since b is bounded on compact sets and b(y) = o(\y\) as y —> oo, for any 6 > 0 
sufficiently small, there exists a constant M > 0 such that |er _1 5(xo + cry)\ 2 < S\y\ 2 + M for any 
y G R d . Thus, 

[ la^bixo + aW^dsK [ (S\W s \ 2 + M)ds <TM+ T6 sup\W S \ 2 

JO Jo s<T 


<TM + TSJ2(( sup Wj) 2 + (inf W° s ) 2 ). 


■ , s<T 
3 =1 


Hence, 


exp ((2p 2 ~P) J W 1 b( x o + crW s )\ 2 ds^ 

d 

< e ( ' 2p2 ~ p ' )MT E exp (t 6{2 p 2 - p) ^((sup W^) 2 + (inf Wj ) 2 )) 


3=1 ^ T 

< e (2p 2 - p )MT( E exp(2T<5(2p 2 -p)|W^| 2 ) ) d/2 , 
where the last inequality follows from Holder’s inequality and the fact that 

sup Wl l = - inf Wj l = \WJf\. 


s<T 


s <T 


Because E 


exp (2TS(2p 2 -p)Wr 


oo. Furthermore, since equation m still holds if one replaces b(x o 
a similar argument yields E[|Zy| p ] < oo. 


< oo if one chooses <5 < (4T 2 (2p 2 —p)) 1 , we obtain E[|Zt| p ] < 

erWs) with b(x 0 + ctW^( s )), 


□ 


Remark 3.3. In general, the conclusion of Lemma \3.2\ is no longer correct if we only suppose that 
b is of linear growth or even Lipschitz. 

Indeed, consider the particular case that d = l,er = 1 and b{x) = x, which is a Lipschitz 
function. It follows from Holder’s inequality that 


E 


exp 


(If 


) E exp( 


exp ( p W s dW s - 


W 2 ds 

= e- pT / 2 E[ex P (| jf W 2 ds'j J E exp (|w 2 - | /' W 2 ds) 


W 2 ds 

r T 
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> 


e-P T / 2 (E[e^/ 4 ]) . 

Furthermore, for any p,T > 0 such that pT > 2 and pT 2 < 1/2, we have E[e pW4r / 4 ] = oo, whereas 


E 


exp 


(ir^) 


< E 

< E 


exp 


SU P \ w s' 2 ' 

v z s<T 

pT 2 pT 


< 


exp (^-(sup W s ) 2 + ^-(mf W 5 ) 2 )) 

(Efe^'^l 2 ]) 2 


Therefore, 


E 


exp 


pT pT 

(p j WsdWs - | J W 2 ds) 


< oo. 


= oo, if pT > 2, pT < 


3.2 Some auxiliary estimates 

We recall that A is a class of functions satisfying (O and W is a d-dimensional Brownian motion. 
The following result plays a crucial role in our argument. 

Lemma 3.4. For any f € A, any p> 1, t > s > 0, 

E[|C(Wt) - CTOH < Cp^Tp, (is) 

V s 

for some constant C p not depending on neither t nor s. On the other hand, if £ is a-Holder 
continuous then 

E[|C(wi) - CTOH < c p (t - s y/ 2 . (17) 

Proof. If £ € A, let (Cn)ngn be an approximating sequence of £ in A. Since £/v —> £ i n ^) oc (K d ) 
and £ and £jv are uniformly exponential bounded, we have 

n\aw t ) - awgm = um eoc^w) - ov(w 5 )h. (is) 

N—>oo 

Next, we will show that 

sup E[|Gv(W t ) - OvTOH < (19) 

JVeN V s 

Indeed, we write 


E[|Civ(m)-Cjv(^)H 

ICiv^ + y) -Cx(a:)| p 


e -\x\ 2 /2s e —\y\ 2 /2(t—s) 


R d J R d 


(27rs) d / 2 (27r(t — s)) d / 2 

<C f [ \C N (x + y)^C N (x)\(e K ^+e K ^)P- 1 j i 
J R d Js. d b 


dxdy 


-\x\ 2 /2s e -\y?/2(t-s) 

(27rs) d / 2 (27r(t — s)) d / 2 


dxdy 
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R d JO 


^j 

<cVt~sf' [ 

J R 


-1 dC, N (x + 6y) 
Vi 


dxi 


w ... . . ,. e -M / 2s e“l y l / 2 (*—«) 

oK(p- l)(\x\+2\y\)f e 


It follows from A(iii) that 


1=1 


r 1 

aCjv(a: + 6»y) 

Jo 


r 

5£7v(a: + %) 

J R d 

dxi 


e -|x| 2 /4s e —|y|/4(t—s) 


(27rs) d / 2 (27r(t — s)) d / 2 
dddxdy. 


dddxdy 


(27rs) d / 2 (27r(t — s)) rf / 2 


0 -|x| 2 /4s 


(27Ts)( d !)/ 2 


da: < Ce^l, 


thus 


E[|Civ(W t )-C J v(W s )n<C 


y/t — S 


Ce K l fl «l 


e -|y| 2 /4(i- s ) 

(2^(t - 


< c 


y/t — S 


From (fT51) and (fl9l) we get (fl6l) . The proof of (H71) is straightforward. □ 

Lemma 3.5. Suppose C,a € A and Qh '■ K d —> K is a-Holder continuous with a € (0,1]. Let M he 
a non-negative constant. Then there exists C > 0 such that for any 0 < t± < £2 < £3 < £4 < T, 


E 


E 


E 


E 


ICaWJ - UWtJWUWtJ - UW^lEeW 


1=1 

4 


ia(w t2 ) - CA(w tl )nc ff (w t4 ) - w 3 )i E eM|lVti 


i=i 

4 


ICff(Wt a ) - CH{W tl )\\U{W ti ) - UW t3 )\ ^e M l^ 


i=l 

4 


ICff(Wi a ) - CffWJIKffW.) - 6KWi 3 )l £ e " |w ‘‘ 


i=l 


< Cy/t\ — to y/t2 — tl 

y/t3~-t2y/ti 

^ C(t4 — t-s) a l 2 y/t- 2 . — t\ 

~ Wi ’ 

< Cyjt\ — to(t2 — t\) a t 2 

y/ts — t2 

<C(t 4 -^) Q/2 (i2-tl)“ /2 . 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


Proof. Let ((a,n)ngn be an approximating sequence of £ a in A. Since (a,n — t (a in L\ oc (R d ) and 
£a and Ca,n are uniformly exponential bounded, we have 


E 


1 C A(Wt a ) - Uw tl )\\U{w u ) - CA(^ 3 )lE eM|WtJ 

1=1 


= lim E 

N—> 00 


4 

ICA.ivWJ - CA,iv(W tl )||CA,iv(W t4 ) - ( A , N (W ti )| Ee M| ^' 

i=l 


(24) 
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Next, we will show that 


sup E 

NGN 


I U, N (W t2 ) - CmtWJIK a,n(W u ) - UAWu)\ E fMlWt * 


< 


Cy/t 4 - t 3 ^/t 2 - ti 


y/ts — 


( 25 ) 


We observe that 
E 


\U, N (W t2 ) - U,N(W tl )\\C A , N (W t4 ) - (A, N (W ti )\J2e MlWtJ 

i= 1 

= / dx dy dz dw\( AtN (z + w) - ( A , N (w)\\( A , N (x + y + z + w) ~ ( A , N (y + z + w)\ 
x {e MM +e M| Z+H +e M| y+ ^| +eM |^^|}^( a .) ft3 _ ta ( y ) fta _ ti ( JB ) fti ( f|> ) 

<c dx dy dz dw\( A , N (z + w) - ( A , N (w)\\( AiN (x + y + z + w) - ( AiN (y + z + w)\ 
J R d JR<i jR<i jR<i 

X 5 c (t 4 -i 3 )( ;c )5c(t3-t2)(y)5 c (t2-ti)( z )ffrf 1 (w). (26) 

Using the mean-value theorem, (l26l) is bounded by 


d( a,n{w + dz) 


dw,. 


d( A , N (y + z + w + 6x) 


dyi 


dx f dy f dz f dw j dd f dS 

l d J R d J R d Js. d Jo Jo 

X 9c(u-t 3 ) ( x)g c (t 3 -t 2 ) {y)9c(t 2 -t 1 ){z)g ctl (w) 

Cy/t 4 — tz\/t 2 — 1\ f dx f dy f dz f dw f dd f dS 
Ju d Jwi d J«. d J«. d Jo Jo 

9c(t i -t 3 ){ x )gc(t 3 -t 2 ){y)9c{t2-t 1 ){z)g c t 1 (w) (27) 


dC a,n(w + dz) 

d( A ,N(y + z + w + Sx) 

dwi 

dyi 


It follows from A(iii) that 


dC, A , n{w + dz) 


/ R d 


dwi 

d( A ,N{y + z + w + dx) 


g 2 ct ^ 


dyi 


(2cnti)( d W 2 

\y\ 2 

g 2 c(* 3—t2) 


( 2 c 7 r(t 3 — t 2 ))*- d 1 ^ 2 


dw < 

dy < Ce K \ z+w+5x \ 


thus (1271) is bounded by 


Cy/t\ — tsy/t2 — t\ 

f 

V^3 ' 

' ^2 

J R d 


Hy/^2 — 

tl 

V^3 

— t2y/ti 


1 

1 

Hy/ti — 

tl 

V^3 

— tzy/ti 



f dz f dw f dd 

d(A,N(w + dz) 

J R d Jm. d Jo 

dwi 


9c(t 4 -t 3 ) (^)ffc(t 2 - tl ) (^)ff ctl (w) 


da; / dzg c{ti -t 3 ){x)g c (t 2 -t 1 ){z) 

JR d 


From (l24l) and (l25l) we get (j20|) . 

The proofs of (17X1) . (l22l) and (l2dl) follow from similar arguments. 


□ 
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In order to prove some key estimations in Lemma 13.71 we need the following inequality. Recall 
that rjh(s) is defined in Jl]). 

Lemma 3.6. Let n be a natural number such that (n — 1 )h < T < nh. Define t/ = ih for 
i = 0,... , n — 1 and t^ = T. Then it holds that 

/■** , / 1 1 1 

Jt’t [VVhi^-Sy/rihis) \Zr]h(s) v r] h (u) — s 


< \\f2 + 2 VT - 


4 ^ 3/2 


E 

2=1 ' 


du 


Proof. We first note that if s > t f { then rjh(s) > s/2. The first integral is estimated as follows 


E 

2=1 ' 


n —1 


du 


ds - 


tf V'thiu) - s-\Jrj h (s) 


< V2hY^ 


pi/ 2 


ds 


ds 


1 


i=2 
n— 1 


< V2h^ 


2=2 


ti / 2 J V— S V$ 

l 

tf/2 Vti - s 


pi/ 2 1 /■* 

—;=ds + 
si* V s 


ds 


< 4a/2. 


The second integral is estimated as follows 

n—1 


E 


i=l “ M “ <-i VVhXpj 
The third integral is estimated as follows 


1 du f ds — -== < y/2h'S~' [ —j= < 2\/T. 

t? M Vm(s) V s 


2-1 /. t ' 1 




ds- 


1 


2-1 




ds 


n—1 


tf Vdhiu) - s i=2 Jt\ vv-s 


2—1 


= 2/*E\A"- /i = 2 E 




d—9 


r T 4T 3 / 2 

< 2 / 5—• 

do 3 


This concludes the proof of the statement. 


□ 


Now by using Lemmas 13.4113.51 and 13.61 we are able to obtain the following estimates on the 
difference of Y t and Y t h , which are defined in Section 13.11 

Lemma 3.7. Let a € (0,1]. Suppose that b £ B(a) defined by (j4]). 

(i) If b is bounded, then for any p > 2, there exists a constant C > 0 such that 

E[| Y r - Yp\ p ] < Chf 12 + Ch pa/2 . 


(ii) If b is of linear growth, then there exists a constant C > 0 such that 

E[| Y t - Yfi | 2 ] < Ch 1/2 + Ch a 
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Proof. Using Minkowski’s inequality, we obtain E[|Yt — Yp | p ] < C{S\(p) + S 2 (p)} where 


Si(p)=E / {(cr 1 b) j (x 0 + aW s ) - (a 1 b) j (x 0 + crTU %(s) )}dW s J 


5 2 (p)=E / {|cr + o-tU s )| 2 - |cr ^(cco + aW Vh ^)\ 2 }ds 

1 Jo 

It follows from Burkholder-Davis-Gundy’s inequality that, 


d r T 


Si(p) < CJ2 / E \(a~ 1 b)j{xo + aW s ) - {a^b)^ + aW Vh{s) )Y 

3=1 Jo 


ds. 


Since b is of linear growth, 


d ph 


J2 E K' 7 ™)j(xo +eWs) - (fT 'fyjixo +crW Vh(s) )\ f 

3 = 1 J ° 

Furthermore, it follows from Proposition 12. II ii) and Lemma HT7T1 that 

[ ^\\i^ 1 b) j (x 0 + crW s )-(a- 1 b) :j {xo+aWr lhis) )\ p 


ds < Ch. 


3 = 1 ' 


ds 


< C 


n/-S , 5 ^ + \ s - ?7n(s)| pa l ds < C(h 1/2 + h pa/2 ). 
VVh{s) ) 


Therefore, Si(p) < C(h 1 / 2 + h pa ). 

Proof of (i). We assume that b is bounded. Using Holder’s inequality, we obtain 


S 2 (p)<E / |cr 1 b(x 0 + crW s )\ 2 - \a 1 b(x 0 + aW Vh ^))\ 2 ds 
l J o 

d 

<CV / E[ \{a~ 1 b)j(xo + aW s )\ 2 — |(<7 _1 6)j(a;o + crU / 7 )h(s) )| 2 H ds. 

3=1 Jo 

Since b is bounded, it holds that for any x,y GM. and j = 1,... ,d, 

|(c r ” 1 6) j (x)| 2 - |(CT- 1 %(y)| 2 P < C{|(cr _1 Mj(:r) - {af L 1 b) j {y)\ p + | {<j- 1 b H ) j {x) - (^ 1 6) j (y)| p } . 

Thus, by dividing the integral into two parts: from 0 to h and from h to T, and applying a similar 
argument as above, we obtain 

S 2 (p) <Ch + C [ ^pMS - ds + Ch pa/2 < C{h l/2 + h pa/2 ). 

Jh VVh{s) 

Thus, E[\Y T - Yp H < Ch 1 / 2 + Ch pa / 2 . 
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Proof of (ii). We assume that b is of linear growth. We observe that 


S 2 (2) = / du dsE {|cr 1 b(x 0 + aW s )\' 2 - \a + ctW i?)i ( s ))| 2 } 

Jo Jo L 

x { \a~ 1 b(xo + aW u )\ 2 - \a~ 1 b(x 0 + crW Vh ( u) )\ 2 } 

Let n be a natural number such that (n — 1 )h < T < nh. Define t= ih for i = 0,..., n — 1 and 
= T. Since b is of linear growth, we have 

S*( 2) 

< Ch + Y^ 


-i ,t h . 


du ■ 


ds 


ds 


ds 


i=l 


xE[{|ct 1 b(x 0 + aW s )\ 2 - \a 1 b{x 0 + aW^^)] 2 } {\cr 1 b(x 0 + aW u )\ 2 - \a 1 b{x Q + aW Vh{u) )\ 2 } 

n — 1 r t^j_i pt h 

<ch+j2 I ' 


du / ds 


i =1 


xE[{|a 1 b(x 0 + aW s )\ 2 - \cr 1 b(x 0 + crW^^)] 2 } {\cr 1 b{x 0 + crW u )\ 2 - \a 1 b(x 0 + <jW Vh{u) )\ 2 } 

Since that for any x, y € R d , 

\a~ 1 b{x) - cr~ 1 b(y)| 2 

d 

= ^2{\(a- 1 b A ) j (x)\ 2 -\(a- 1 b A ) j (y)\ 2 

.7=1 

+ {(a-'b^jix) + {(p- 1 b H )j(y)}{(<J- 1 bH)j{x) - {(v~ 1 b H )j(y)} 

+ 2{a~ 1 b H ) j {x){{(j- 1 b A ) j {x) - (a~ 1 b A )j(y)} + 2{a- 1 b A ) j (y){(a~ 1 b H ) j (x) - (cr _1 ^)j(?/)}}. 
by using Lemma [3751 with t A = rjh(s), t 2 = s, = r]h(u) and 1 4 = it, we obtain 


S 2 { 2 ) < Ch + Cj2 


— 1 


2=1 


du ds ■ 
It’y Jt 5 * 


Vvh{u) - SyJy h {s) \fryjs) \/vh(u) - s 


Therefore from Lemma 13.61 we have 

S 2 (2) < C(h + h* + % +h a ). 
Thus, E[|y r - Y$ I 2 ] < Ch V 2 + Ch a . 


□ 


3.3 Proof of Theorem 12.51 

We first recall that Z t , , Y t and Y t h were defined in Section [3TT1 It follows from Lemma [3771 that 
Yj f converges in probability to Yt as h —> 0. Thus Z ^ also converges in probability to Zt as h —> 0. 
Moreover E[Zy] = E[Zt] = 1 for all h > 0. Therefore, it follows from Proposition 4.12 in ITS] that 

lim E[|Z£ - Z T |] = 0. (28) 
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On the other hand, since / is bounded, it follows from m that 

\E[f(X)~f(X h )]\<CE{\Z£-Z T \\. 

This estimate together with C51) implies the desired result. □ 

3.4 Proof of Theorem 12.61 

It is clear that \e x — e v \ < ( e x + e v )\x — y\. This estimate and Holder’s inequality imply that 
|E {f(X) - f(X h )\ | is bounded by 

E[\f(x 0 + oW){Z T + Z$){Y t - Y£)\] 

< \\f(x 0 + aW)(Z T + Z£)\\ 2 \\Y T - Y £|| 2 

/ \ 2/r / \ (r—2)/r 

<{E[\f(x 0 + aW)\ r ]) [E[\Z t + Z^\ 2r/< ~ r ~ 2) ]j \\Y t -Y£\\ 2 . 

Thanks to the integrability condition of / and Lemma 13.21 

(e[|/(x 0 + aW)\ r ]y / ' (e[\Z t + Z^\ 2r ^ r - 2 ^y r ~ 2)/r < C < 00 . 

This together with Lemma 13.71 implies the desired result. □ 


3.5 Proof of Corollary 12.71 

Since g is /3-Holder continuous, it holds that E \\g (max 0 < s <T |*o +crW s |) | r ] < 00 for any r > 2. 
Thanks to Theorem 12.61 it remains to estimate 


E 


9 


max 

0<s<T s 



I max \X 

\0<s<T 



Since g is /3-Holder continuous and b is sub-linear growth, we have 


E 


q I max \X, 

a '0<s<T 1 s 


-E 


q max 

'0<s<T 


\X h 


Vh(s) 






r -1 

< CE 

max \Xg | — max IJ'C 1 

0<s<T sl 0<s<T Vh{s)l 


< CE 

max \X% - X% ( ) \ 

0<s<T s rlh(s> 


< CE 


o max r {|6(^ (s) )|( S - r] h (s)) + \a{W s 


W Vh ( s) )\} 


< Ch p + CE 


max I W s 

0 <s<T 


W Vhi s) 1} 


By modulus continuity of Brownian motion (e.g. Lemma 4.4 in E3). we have 


E 


max I Wt 

0 <s<t<T,\t-s\<h 



< Ch log (1 /h ), 


( 29 ) 
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Thus from (1^1) . we obtain 


E 


9 


max IJC 

0 <s<T 



( max 

0 <s<T 


\X h 


Vh(s) 


<c[h^+ {h\og{l/h)f/ 2 Y 


which implies the proof of the statement. 


□ 


3.6 Proof of Corollary 12.81 

We first note that if b is bounded, then it holds from Theorem 2.1 in [26] (see also Corollary 3.2 
in [33]) that there exists a density function p of for t £ (0,T] and it satisfies the following 
Gaussian upper bound, i.e., 


p h t (x) < C- 


.Mor 


t d / 2 


for some positive constants C and c. 


Now we prove that E 


it follows from Jensen’s inequality that for any r > 2, 


/ ^ Jq g(xo + W s )ds) is finite for any r > 2. Since |g r (a;)| < Ke K ^, 


E 


/ / g(x o + w s )ds 


r T d 


<C + C E[\g(x + W s )\ r ]ds<C + C \\E[e rKW °]ds. 


Since x 2 /(As) + K 2 r 2 s > Krx , we have 

r T 


E 


f l J g{xo + w s )ds 


<c + c 


„dKrs 


Thanks to Theorem 12.61 it remains to prove that 


E 

/ ( [ T g(X^)ds) 

-E 

/ ( r^ds) 


l V° J\ 


\Jo J J 


ds < oo. 


< Ch p/2 . 


Since / is a Lipschitz continuous function, we have 


E 


/( g(X^dsj^-E^f^ g(X^ (s) )dsj j <Cj^ E[g(X!)-g(X* kW ) 
If s 6 (0, h], then by using the Gaussian upper bound for Pg(x), we have 


ds 


E 


g{X d )-g{ Xl (s) ) ds< E[\g(X^\]ds + \g(x 0 )\h 
J Jo 

^ I X —XQ | 

< C [ ds j \g{x )\-—^— + \g{x 0 )\h 
JO J R d s ' 
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< Ch. 


On the other hand, for s G [h,T\, using the Gaussian upper bound for P^ h ^ and following the 
proof of Lemma [3~~il (see also Lemma 3.5 in [S3); we have 

E[| 9(X!) - g(X v h h(s) ) |] < + Ch^\ 

VVhis) 

Therefore, we conclude the proof of the statement. □ 

3.7 Proof of Corollary 12.111 

Recall that r|^, r^’ h , Z , Z h were defined at the beginning of Section [3.11 It suffices to proof 
the statement for the case that g is positive. It follows from (fl2l) that E[g(AT)l( rD >T)] — 
E[#(X£;)1( t £ >T )] = Ei + E 2 where 


Ei = E[g(a;o + aWr){Z T — Z^)l^w,h >T ^], 

E 2 = E[g(a;o + aW T )Z T {l^w >T) - 1 ( r w,h >T) )]. 

It follows from the proof of Theorem 12.61 that 

|^i| < E[|«?(z 0 + <tWt)(Z t - Z %)|] < C7i* a *. 

Applying Holder’s inequality, we have 


\E 2 \ < \\Z T \\ q \\g{x 0 + CrWr)(l( r w >T) - l( T w,h >T) )||p, 


where q is the conjugate of p. Thanks to Lemma 13.21 and the fact ,h > , we have 


) 1 (r^ h >T) 


-E 


g p (x 0 + crWr)l (T w> T) ) 


\E 2 \<C p (E[gP(x 0 + aW T )l 
It follows from Theorem 2.4 in [TO] that there exists a constant K(T) such that 

/K(T)\W’\UV / \± 

|£2| - C H 1A e 4 I h ■ 

Combining this estimate with ( 13 ( 71 ) completes the proof. 


i /p 


(30) 


□ 


3.8 Proof of Theorem 12.121 


The proof is the based on the perturbation or Levi’s parametrix method (see 0) for the density 
functions pt(xg,y) and Pt(xo,y). It is known that when the drift coefficient b is bounded and 
Holder continuous, it holds that for any t € (0,T] 


Pt(xo,y) = gta{x 0 ,y) + 



v*g«-,)a{z, y),H z ))pt{x o, z) 
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( 31 ) 


= gta{x 0 ,y) + / E [(X7 x g( t -s)a{X s ,y),b(X s ))] ds, 

Jo 

where a = a a* and g ta (xo,y) is called the parametrix (see, Chapter 1 in j7]). 

We first consider a representation which is similar to (1311) for the density functions pt{xo,y) 
and Pt (xo, y) under the assumption that the drift coefficient is bounded measurable. Recently, 
Makhlouf [BUI Theorem 3.1] proved that the representation (13T1) also holds for a Brownian motion 
with random drift b = ( bt)o<t<T under the suitable growth condition. For the convenience of the 
reader, we give a proof below. 

Proposition 3.8 (Makhlouf [3D] ) . Let W = {Wt)o<t<T be a d-dimensional (fFt)o<t<T-Brownian 
motion. Suppose that the stochastic process b = ( bt)o<t<T is adapted to ( jFt)o<t<T o,nd there exists 
a constant K > 0 such that sup 0<t<T |6 t | < K almost surely. Then, for any t £ (0,T] and x £ R d , 
the stochastic process Y t := x + j^b s ds + aWt admits a density jt(x, •) with respect to Lebesgue 
measure and for any y £ R d , 


7 t(x,y) 


9ta(x,y)+ / E [(V x g {t _ s)a (Y s ,y),b s )] ds. 
Jo 


(32) 


Proof. It suffices to prove that for any infinitely differentiable functions / : R d —> R with compact 
support 

d yf(y) [ dsE 

Jo 

It is well-known that the function u(s, x) := E [/ (x + oWt~ s )\ is a solution to the following partial 
differential equation: 


— s)o(W; y\ ^s)] ■ 


E [fO / t)]= [ f(y) 9 ta{x,y)dy+ [ 

J R d jR d 


„ , ,1 \ ^ cl 2 , . 

9 s u(s,x) + - 2^ a id -—— : u(s,x) 


i,j =1 


dxidxj 


u(t , x) 


0, (s, x) £ [0, t) x R d , 
/(x), x £ R d . 


(33) 


Hence we have 


E [/ (x + crWt)\ = u(0, x), 

E[f(Y t )\=E[u(t,Y t )]. 


By using Ito’s formula and (1331) . it holds that for any e £ (0, T), 

u(t-£,Y t - e ) = u(0,x)+ f (\7 x u(s,Y s ),b s )ds + ^2 I a itJ1 J-u{s,Y s )dWf. 

Jo 1 Jo 


i,j =1 ' 


Since for any i = 1, ..., d, s G [0, T) and x G M d , 

) 9 f d 

-u(s,x) = — E[f (x + aWt-s)} = J f{y)-Q^g(t-a)a{x,y)dy 


_d_ 

dxi 


(34) 

(35) 


(36) 


(37) 
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and 


9 ( \ 
—u{s,x) 

UXn 


< 


cWfU 


for some constant C > 0, the stochastic integral in (l36l) is martingale. By taking the expectation 
and Fubini’s theorem, we have from (13-H) 

E[u(t - e,lt_ e )] = E [/ (x + aW t )} + [ E [{V x u(s,Y s ),b s )]ds. 

J o 

Taking e —> 0 and using the dominated convergence theorem, we have from (1351) and (1371) . 
E[/(it)] = lim E‘[u(t-e,Y t - s )]=E[f(x + o-W t )]+ [ E [(S7 x u(s, Y s ), 6 S )] ds 

e ^o+ J 0 

= [ f{y)dta{x,y)dy+ [ f(y) f E [(V x g( t - s)a (Y s ,y),b s )] dsdy. 

JK. d J«. d JO 

This concludes the proof. 

Proof of Theorem \2.1°A Using Proposition 13.81 we have 

p T (xo,y) = g T a{x 0 ,y) + / E [(V x g^ T -s)a(Xs,y),b{X s ))] ds, 

Jo 


□ 


Pr(xo,y) = gTa(x 0 ,y) + E CV x g( T - s )a{X s , y), b{X Vh{s) ) 


ds. 


Moreover, from Lemma o we have 


p T (x Q ,y) - PT(x 0 ,y) = [ E [{V x g^ T _ s ) a (x 0 + crW s ,y ), Z T b(x 0 + aW s ) - Z^b(x 0 + a w ))] ds. 
Jo 

By using Jensen’s inequality and Schwarz’s inequality, there exists c > 0 such that 

\pt(xo,v) -Pt( x o,v)\ < f E [|V x 3 (T _ s)a (a:o + W s ,y)\ \Z T b(x 0 +aW s ) - Ztfb(x 0 + CTW„ h(g ))|] ds 

Jo 


< C 


< C 


VT- 


=E [5c(t- s )(zo + W 8 ,y) \Z T b(x 0 + aW s ) - Z!fb(x 0 + |] ds 




- (^4 S + B s )ds , 


where 


A s := E [g c (T-s)( x o + aW s ,y) \Z T - Zj*\] 

B s := E [g c (T- s ){x 0 + aW s ,y ) \Z T \ \b{x 0 + aW s ) - b(x 0 + <xW Vh ^)\ 
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Note that for any q > 1, there exist C q and c q such that 


E[\g c{ T- s )(xo + aW s ,y)\^} 1/q < C, 


T- 


d(q-l) 

2 q 


9c q T(x 0 ,y). 


(38) 


We first consider the upper bounded for A s . By using \e x — e v \ < (e x + e v )\x — y\ and Holder’s 
inequality, for any p, q > 1 with l/p+l/q = 1 and r > 1, we have 


< E [\g c{T - s) (xo + aW s ,y)\«] 1/q E \\Z T + Z^\ p \Y t — Yr 


rh\P 
1 T 


1/p 


<c q 


d(g—1) 

T \ 2q 


,T- S/ 
Therefore, it holds that 

rT 


gc q Hxo,y)E\\Z T + Z^r/^ 


(r-l)/(pr) 


E 


i Y T - Y r r 


■1 1 /(pr) 


/0 


< C p , r C q g CqT (xo,y) J o VT — ( T _ g 


d(q-i ) 

T \ 2q 


dsE 


I Ft-KtT 


1 /(.pr) 


• (39) 


By choosing p > d that is q = p/(p — 1) < dj (d— 1), from Lemma 1X71 we have 


\Yt-Y t T 


A s ds < C p , r C q g c T (x 0 ,y)T 1/2 E 

V J- — s 

< C p ^ r g Cq T(.x 0 ,y)T l / 2 {h a > 2 + h 1 '^}. 


1 /(pr) 


(40) 


Now we give an upper bound for B s . By using Holder’s inequality, for any p > d and 
q = p/(p — 1) and r > 1, we have 


B s < E [\g ciT -s)(xo + aW s ,y)\ r >] 1/q E [|Dl v ' ^"" 'e [| b(x 0 + aW s ) - b(x 0 + aW Vh{s) )\ pr ] 


(r-l )/(pr) 


1 1 /{pr) 


— ^PiQi r 


d(q-l) 

T \ 2q 


K T — s 

By Lemma 13.41 for any s > h, 


g Cq T(x 0 , y )E [|fo(ar 0 + crW s ) - b(x 0 + aW Vh ( s )) 


I pr -i l/(pr) 


111 \ 1,1/( 2 pr) 

E [|6(*o + vW.) b{x o + *W vkW )H ' (pr) < C{h^ + Vh{s)1/{2pr) }- 


Since b is bounded, we have 

rT 


l VT — s B ° ds ^ c p, q xgc qT (x< h y)J o VT —( T s 

+ C p ^ r g CqT {xo,y)h a/2 —l 


d(q-l ) 

T \ 2q 


ds 


d{q- 1 ) 

T \ 2q 


-s \T~S 


ds 


C p ^ r g CqT {x 0 ,y)h}/^ 


VT~s \T — s 


d(q-l) 

T \ 2q 


Vh 


Js) (41) 
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By the assumption h € (0,T/2), it holds that 


VT— \T- 


d(q-l) 

T \ 2q 


1 i d(q-l) fl 

ds < 2 2 + ^5— 


T l/2 ■ 


Since q < d/(d— 1), we have 


VT~s \T — s 


d(g-l) 

1 \ 2q 


ds = 


IQ , . g —dg+d 

(T -h) 2 « 


q — dq + d 


and 


d(g-i) 
1 \ 2g 


Jh VT- s \T-sJ %(s) 1/(2p ' 
where B(x,y) is the beta function. Therefore, we obtain 


i —-ds <(T- h) 1 - 11 ^ b(- + 1} , —") 

i/( 2 P T-) - > ^2 2q 2prJ 


Vtt s 

which concludes the proof. 


£ ^f^= s B ° ds < C p , q ,r9cT{x 0 , y ) + ^ /2 + , 


(42) 


□ 


3.9 Proof of Theorem 12.161 

In the same way as in subsection ld.il we have the following Lemma. 

Lemma 3.9. Let X be a solution of the SDE (j8]) and let X h be the Euler-Maruyama approximation 
of X defined in © with h > 0. Suppose that b is a measurable function with sub-linear growth. 
Let U = ( Ut)o<t<T be the unique solution of the equation Ut = x o + crWt + L®(U) and 

Z t := e*, Y t := f b(U s )d.W s - \ f b 2 {U s )ds 
Jo ^ Jo 

Z? := Y t h := £ b(U Vh{s) )dW s - ± j* b 2 (U Vh{s) )ds. 


Then it holds that 


E[f(X)} - E[f(X h )] = E [f(U)(Z T - Z%)\ 

for all measurable functions f : C[0, T] —> K. provided that the above expectations are integrable. 
Proof. We define new measures Q and Q h as 

^=exp J q cr~ 1 b(X s )dW s \(r- 1 b{X a )\ 2 ds'), 
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dp 


eXP ( " Jo a ~^ X l(°) )dW ° - \ Jj \^ X t(s))\ 2 ds) ■ 


Since b is of sub-linear growth and the fact that 0 < L^(X) < \a\ sup 0<s<4 \W S \ 1 by following the 
proof of Lemma HO we can show that Q and Q h are probability measures. Furthermore, it follows 
from the Girsanov theorem that the processes B = (B t )o<t<T and B h = {B^) 0 < t < T defined by 


Bt = W t + / a^biX^ds, B* = W t + a-^X^Ms, 0 < t < T, 

Jo Jo 

are Brownian motions with respect to Q and Q h respectively. Note that X s = Xq + crB s + L° S (X) 
and X'f = xo + <xB^ + L° s (X h ). Therefore, 


E[.f(X)} = E d f(X) 


dP 


- e q 
= Ea 


1 


/(X)exp(J( a~ 1 b(X s )dB s -^ \a~H{X a )\ 2 da^ 
f(xo + crB + L°(X)) exp ( / a~ 1 b(x 0 + uB s + L° s (X))dB s - ^ [ \a~ 1 b(xo + crB s + L° s (X))\ 2 ds 


Since (X , B) |q = (U , TT)|p, the above term equals to 

r T 


E 


1 


f(xo + crW + L U (U)) exp ( / a 1 b(xo + aW s + L^(U))dW s - ^ J \a 1 b(x 0 + oW s + L° s (U))\ 2 ds 


= E 


f{U) exp ( / a~ 1 b(U s )dW s — ^ ^(U^ds) 


= E [f{U)Z T ]. 

Repeating the previous argument leads to E [f(X h )] = E[f(U)Zj], which concludes the statement. 

□ 

In the same way as Lemma 13.21 we have the following estimate of the moments of Z and Z h . 

Lemma 3.10. Suppose that b is of sub-linear growth. Then for any p > 0, 

E[\Z T \ P + \Zlf\P\ < C < oo, 

for some constant C which is not depend on h. 

Finally, we introduce the following auxiliary estimate. 

Lemma 3.11. Let U as in Lemma \3.tA Suppose that f is a-Holder continuous with a € (0,1], 
then for any t > s > 0, 

E[\{(u t )-aus)\ p ]<c p (t-sy a / 2 . 
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Proof. By Holder continuity of we have 

E[| <Z(Ut) - C(U s )n] < E [|Ut - u s \ ap ] < CE[\W t - W S H + CE[\L° t (JJ) - L°(U) H- 

Hence it is sufficient to prove that 

E[| L° t (U) - L° s (U)\ pa ] < C p (t - s) pa/2 . 

Using Lemma r2.1dl we have 

L° S (U) < L° t {U) < L° S (U) + sup max (0, —a(W u — W s )). 

s<n<t 

Therefore, \L°(U) — L° S (JJ)\ < |tr| sup s<u<t | W u — W s \. Hence applying Burkholder-Davis-Gundy’s 
inequality, we have 

E[| L° t (U) - L° s {U)\ pa } < CE[ sup | W u - W s \ pa ] < C(t - s) pa / 2 . 

s<u<t 


This concludes the proof. 


□ 


Proof of Theorem \2.16\ We first recall that Z t . Zf , Y t and Y t h are defined in Lemma [3.91 Using 
Lemmas 13.91 and the elementary estimate \e x — e y \ < ( e x + e v )\x — y |, we have 


|E lf(X)} - E[f(X h )}\ < E | f(U)(Z T + Z*)(Y t 



Thanks to Lemma T3.101 and the Holder’s inequality, for some r > 2, we have 
|E[/(X)] - E[f(X h )]| < CE[\f(U)\ r ] 2/r \\Y T - Y £|| 2 . 


By a similar argument as the proof of Lemma 13.71 we can show that 


\\Y T -Y£\\ 2 <Ch a/2 , 


which concludes the proof. □ 

Remark 3.12. The conclusion of Theorem \2. 16\ still holds if we relax the condition f bounded to 
E[|/(C/)| r ] < oo for some r > 2. 
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